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The  interaction  between  an  intense  radiation  field  and  atoms  that  undergo  collisions  in  an  atomic 
vapor  is  studied.  With  the  assumption  that  velocity-changing  collisions  are  state  Hide  undent  and 
weak,  w<j  obtain  an  analytic  expression  for  (he  time-dependent  density-matrix  elements  in  solving  a 
quantum-mechanical  transport  equation  using  a  dressed-state  approach  Our  analytical  calculation 
is  used  to  predict  the  form  of  the  signal  in  an  extended-pulse  photon  echo  The  results  show  good 
agreement  with  an  experiment  carried  out  by  Yodh  et  at.  [Phys.  Rev.  Lett  53,  659  M984  )  It  is 
found  that  the  decay  rate  of  the  echo  signal  decreases  rapidly  with  increasing  duration  of  the  second 
pulse  in  the  pulse  excitation  sequence.  A  comparison  between  theory  and  experiment  is  presented 
and  a  physical  interpretation  of  the  results  is  given. 


1.  INTRODUCTION 

The  interaction  between  radiation  and  matter  involv¬ 
ing  various  relaxation  processes  in  gases  and  solids  has 
been  studied  extensively  both  theoretically  and  experi¬ 
mentally  using  optical  coherent  transient  techniques. 
Transient  experiments,  such  as  the  photon  echo,  the 
stimulated  echo,  and  optical  nutation  can  be  particularly 
useful  in  investigating  atomic  collisional  processes.1  111 
Our  present  work  relates  directly  to  an  experiment  car¬ 
ried  out  by  Yodh  et  at. 11  using  a  variation  of  the  conven¬ 
tional  photon-echo  technique.  The  general  methods 
developed  here  can  be  further  used  for  studying  relaxa¬ 
tion  phenomena  in  other  transient  experiments. 

In  the  experiment  of  Yodh  et  at.,  a  two-pulse  excita¬ 
tion  scheme  was  used  to  produce  the  echo.  In  contrast  to 
the  excitation  scheme  of  a  conventional  photon  echo,  the 
duration  of  the  second  pulse  in  this  extended-pulse  pho¬ 
ton  echo  (EPPE)  was  comparable  to  the  atomic  lifetime. 
The  collisional  decay  rate  in  this  experiment  is  found  to 
decrease  rapidly  when  the  duration  of  the  second  pulse 
increases,  an  effect  which  cannot  be  correctly  accounted 
for  by  the  conventional  optical  Bloch  theory.11 

It  is  the  purpose  of  this  paper  to  analyze  the  experi¬ 
mental  result.  The  physical  system  under  consideration 
consists  of  an  ensemble  of  two-level  atoms  (optically  ac¬ 
tive)  that  interact  with  an  external  radiation  field  and  col¬ 
lide  with  perturber  atoms.  Over  the  last  20  years,  pro¬ 
gress  has  been  made  in  two  parallel,  but  independent  ap¬ 
proaches  to  modeling  physical  systems  involving  stochas¬ 
tic  processes  such  as  velocity-changing  collisions.  The 
first  one,  in  which  the  atomic  velocity  appears  in  the 
Bloch  equations1’  as  a  stochastic  variable,  is  somewhat 
more  general.  After  solving  the  Bloch  equations  for  the 
density-matrix  elements  as  a  function  of  velocity  and 
time,  a  statistical  average  over  the  whole  collision  history 
must  be  made.  With  this  approach,  it  is  possible  to 
choose  various  stochastic  models  characterized  by  corre¬ 
lation  functions  describing  different  types  of  fluctuations. 
However,  when  solving  transient  problems,  such  statisti- 
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cal  averaging  procedure  can  he  extremely  difficult.'  A 
discussion  of  this  method  is  not  the  subject  of  the  work. 
The  second  approach  rests  on  a  critical  assumption  that 
the  stochastic  fluctuation  arising  from  velocity-changing 
collisions  m  a  vapor  can  be  described  by  an  impact,  or 
Markovian  approximation.14  namely,  the  atom-perturber 
interaction  occurs  on  a  time  scale  that  is  instantaneous 
with  respect  to  all  other  time  scales  in  the  problem  except 
for  the  optical  period.  The  velocity  of  an  atom  depends, 
at  most,  on  its  value  immediately  before  a  collision.  In 
this  spirit,  the  statistical  averaging  procedure  can  be  for¬ 
mulated  in  a  manner  that  leads  to  a  transport  equation  in 
which  the  stochastic  process  of  velocity-changing  col¬ 
lisions  is  characterized  by  both  an  additional  decay  term 
and  an  integral  source  term.  All  variables  appearing  in 
the  equations  are  well-defined  quantum-mechanical  quan¬ 
tities;  consequently,  the  equation  has  been  termed  a 
“quantum-mechanical  transport  equation"  (QMTE).15  20 
As  in  the  theory  without  collisions,  the  echo  signal  can  be 
obtained  by  solving  the  QMTE  in  each  time  interval  of 
the  excitation  scheme.  Essentially,  the  QMTE  can  be  de¬ 
rived  in  two  representations,  using  either  a  bare-state  pic¬ 
ture  (BSP),21  24  or  a  dressed-state  picture  (DSP).25'26  In 
the  BSP,  the  atom-field  interactions  are  represented  in 
terms  of  a  basis  using  atomic  eigenstates.  In  the  DSP,  a 
new  basis  is  formed  using  atom-field  eigenstates.  In  gen¬ 
eral  (for  arbitrary  field  intensity  and  collision  kernels),  the 
QMTE  in  both  representations  remains  in  a  coupled 
differentio-integral  form  which  is  very  difficult  to  solve. 
However,  for  certain  types  of  physical  systems  encoun¬ 
tered  in  many  experiments,  further  approximations  and 
assumptions  with  respect  to  the  field  intensity,  duration 
of  excitation  pulses,  and  collision  models  can  be  made. 
As  a  result,  the  quantum-mechanical  transport  equation 
is  simplified  and  some  analytic  solutions  can  be  obtained. 

Previous  work  has  shown  that  echo  problems  can  be 
solved  by  applying  the  QMTE  in  the  following  limiting 
cases. 

(1)  The  role  of  velocity-changing  collisions  on  photon 
echoes  has  been  studied  by  a  number  of  authors'  for 
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excitation  pulses  which  are  of  short  duration  [hereafter 
referred  to  as  conventional  photon  echoes  (CPE)].  In  this 
approach,  all  relaxation  effects  can  be  neglected  during 
the  atom-field  interactions.  In  field-free  time  intervals, 
the  QMTE  can  be  solved  for  a  specific  choice  of  collision 
kernel. 

(2)  In  contrast  to  the  conventional  photon  echo,  the 
second  pulse  is  extended  to  be  comparable  to  the  atomic 
lifetime  in  an  extended-pulse  photon  echo  (EPPE).  Re¬ 
laxation  effects  (including  velocity-changing  collisions) 
can  no  longer  be  neglected  in  the  second  pulse  interval. 
To  obtain  the  echo  amplitude,  one  must  solve  the  coupled 
differentio-integral  form  of  QMTE  in  the  second  pulse  re¬ 
gion,  for  which  no  simple  analytical  solution  is  available. 
A  formal  method  for  solving  the  QMTE,  expressed  as  a 
series  in  powers  of  the  excitation  field  amplitude,  has 
been  given.19  The  method  is  useful  practically  when  the 
order  of  iteration  is  low,  implying  that  the  field  strength 
(in  frequency  units)  is  small  compared  with  the  inverse 
pulse  duration  of  the  excitation  field.  We  applied  the  per¬ 
turbation  solution  of  the  QMTE  to  the  extended-pulse 
photon  echo  in  a  comparison  paper.11  Although  the  per¬ 
turbative  picture  does  not  usually  mirror  the  experimen¬ 
tal  situation,  it  can  help  one  to  view  the  underlying  phys¬ 
ics  of  the  collisional  modification  of  the  echo  signal  as  a 
function  of  the  second  pulse’s  duration. 

(3)  Kryszewski  and  Nienhuis’2  have  discussed  a 
method  for  solving  the  QMTE  in  the  study  of  light- 
induced  drift  for  a  somewhat  more  general  collision  ker¬ 
nel  proposed  by  Keilson  and  Storer"  that  can  model  both 
“weak”  and  “strong”  collisions.  The  idea  of  their 
method  is  to  expand  the  density-matrix  elements  in  terms 
of  eigenfunctions'4  of  the  Keilson-Storer  kernel  and  to 
solve  the  resulting  equations  for  the  expansion 
coefficients.  Although  there  is  no  restriction  on  field  in¬ 
tensities  for  employing  the  method,  the  solution  remains 
in  the  form  of  a  recurrence  equation  in  Laplace  transform 
space  that  can  be  solved  only  numerically.  With  a 
specific  initial  condition,  the  solution  is  in  the  form  of  a 
continued  fraction  which  also  requires  a  nontrivial  nu¬ 
merical  computation. 

In  this  paper,  we  discuss  a  new  method  for  solving  the 
QMTE  analytically  for  the  strong-field  limit  and  apply 
the  result  to  the  extended-pulse  photon  echo,  In  Sec.  II, 
the  approximations  and  assumptions  of  the  theory  are 
discussed.  The  relationship  between  dressed-  and  bare- 
state  representations  is  established,  and  the  QMTE  is 
written  in  both  the  bare-state  (BSP)  and  dressed-state 
(DSP)  pictures.  In  Sec.  Ill,  a  solution  of  the  QMTE  in 
the  dressed-state  picture  in  the  strong-field  regime  is  ob¬ 
tained  and  the  solution  is  used  to  obtain  an  expression  for 
extended-pulse  photon  echo  amplitude  in  Sec  IV.  A 
comparison  between  the  strong-field  result  and  the  exper¬ 
iment  is  given  in  Sec.  V. 

II.  APPROXIMATIONS  AND  FORMALISM 
A.  Basic  theory  and  equations  in  the  BSP 

The  physical  system  under  consideration  consists  of  an 
ensemble  of  two-level  atoms  immersed  in  a  vapor  consist¬ 
ing  of  perturber  atoms.  The  active  atoms  interact  with 


the  laser  field  while  undergoing  collisions  with  perturber 
atoms  (see  Fig.  1).  We  make  the  following  approxima¬ 
tions  and  assumptions  with  regard  to  characteristics  of 
the  laser  light,  active  atoms,  perturber  atoms,  and  col¬ 
lisional  processes. 

(1)  The  laser  light  is  taken  to  be  of  the  form 

El{2,t)  =  X£0it)cos(k2 -ill)  (2.1) 

polarized  in  the  %  direction  and  propagating  in  the  2 
direction  with  propagation  vector  k  =  (Il/c)Z.  The  en¬ 
velope  function  E0(t)  is  slowly  varying  compared  with 
cosfir. 

(2)  The  applied  field  frequency  11  is  chosen  to  be  reso¬ 
nant  with  the  atomic  transition  frequency.  In  the  pres¬ 
ence  of  spontaneous  emission,  but  neglecting  collisions, 
atomic-state  density-matrix  elements  evolve  as  (in  a  field- 
interaction  representation) 

pn(v,f  )=  -y,pn(v,i  )  +  y: ) 

_l  V[pi:*vT  )~P;|(v,f )]  .  (2.2a) 

p22(  y,t  )  =  —  )  +  /'Y[pii(  y,f )  — p,,(  v,t )]  ,  (2.2b) 

Pi2(v,/ )=  -(yl2-ifcv)-i:g[p,|(v,t  )-p22(v,r  I]  ,  (2.2c) 

p2,(v,r  )=  -( j'j,  -t  ikv)p2i(v,i  )  +  i\'[fju(vj  )  -p,,< »-,/)]  , 

(2. 2d) 

where  \  is  a  Rabi  frequency,  given  by 


y,  (/=  1,2)  is  the  spontaneous  emission  rate  of  level  i,  y,  , 
is  the  spontaneous  relaxation  rate  from  level  2  to  level  1, 
Yi2  —  y2i=(yi+y2)/2,  and  p  1 2  is  a  dipole  moment  ma¬ 
trix  element. 

(3)  With  regard  to  collisions,  the  following  assumptions 
are  made. 

(a)  Collisions  are  not  sufficiently  energetic  to  induce 
transitions  between  states  1  and  2. 

(b)  The  duration  of  a  collision  r(.  is  assumed  to  be 
smaller  than  any  (i)  relaxation  times  (i.e. ,  }'l2rt.  «1, 


FIG.  1.  Atomic  vapor  consisting  of  two-level  active  atoms 
and  perturber  atoms.  The  active  atoms  interact  with  laser  light 
having  an  oscillation  frequency  11  and  Rabi  frequency  y.  and 
collide  with  perturber  atoms. 
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r rc  « 1 ,  r=collision  rate);  (ii)  the  time  that  the  field 
takes  to  induce  atomic  transitions,  i.e.,  yrc  «  1;  (iii)  the 
inhomogeneous  dephasing  time  experienced  by  atomic  di¬ 
poles,  i.e.,  kurc  «  1,  where  u  is  the  most  probable  speed 
of  the  active  atoms.  These  assumptions  collectively  are 
known  as  the  impact  approximation  (Markovian 
approximation)  —  the  atom-field  and  collisional  interac¬ 
tions  contribute  independently  to  p  . 

(c)  The  density  of  active  atoms  is  low  enough  to  neglect 
active-atom-active-atom  collisions,  so  that  the  only  col¬ 
lisions  that  need  be  considered  are  active- 
atom-perturber-atom  collisions. 

(d)  The  perturber  reservoir  is  sufficiently  large  to  be 
unaffected  by  collisions.  This  assumption  is  implied  in 
(c). 

Thus,  the  collisional  contribution  to  the  time  rate  of 
change  of  density-matrix  elements  can  be  added  to  Eqs. 
(2.2)  to  arrive  at  the  QMTE 

Pn<v,f)  =  -ir->-ylipu(v,t  >  +  r2  ,iP22(v./) 

-'V[Pl2(v,f  )~p21(v,f  )] 

+  f  ui(v'  —  v)p||(v',f  )dv'  ,  (2.4a) 

p22(v,r  )  =  -  (  r  T  y  2  )p 2 2 f  v, 1  +' V[pi2<v.f )  p2 1 ( v,  t )  ] 
+  f  w(  v'  —  v)p22(  v )dv‘  ,  (2.4b) 

pl2(v,r)=-(T  +Tph  + yl2-wv)p12(v,r ) 

-'TtPii(v,/)-p22*v,f )] 

+  f  w(v'  —  v)p12(v',r  Wv'  ,  (2.4c) 

p2l(  v,t  )=  —  (  r  +  Tph  +  y,2  +ik  v)p2l(  v,t ) 

+  '¥[Pn<vT>-p22(v,/)] 

4-  J*  uj(  v'  —  v)p21(  v',r  )dv‘  .  (2.4d) 

In  writing  Eq.  (2.4)  several  assumptions  regarding  the  na¬ 
ture  of  the  collisional  interaction  have  been  incorporated. 
The  velocity-changing  collisions  are  assumed  to  affect 
level  populations  and  coherence  in  the  same  manner. 
These  collisions  are  characterized  by  collision  kernel 
w{v'~v)  and  rate  T.  The  parameter  Tph  in  Eq.  (2.4) 
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represents  a  phase-interrupting  collision  rate  which 
affects  atomic  coherence  only.  The  detailed  description 
of  the  model  and  the  validity  conditions  can  be  found  in 
our  companion  paper  for  the  weak-field  limit. 

The  collision  kernel  to  be  used  in  this  work  is  a 
“difference"  kernel,  given  by 

Uiv'  -vl-  r _  e  >'v  (2.5) 

a  v7  TT 

that  has  been  proven  to  be  a  good  model  for  characteriz¬ 
ing  weak  collision  processes. |l)  The  quantity  a  is  roughly 
the  rms  change  in  velocity  per  collision  which  is  assumed 
to  be  much  less  than  the  most  probable  speed  u. 

Notice  that  in  Eq.  (2.4)  the  Rabi  frequency  serves  as  a 
coupling  constant  which  complicates  the  calculation 
For  short-duration  pulses  or  weak-field  strengths  such 
that  the  pulse  area  lyT^  )  is  less  than  unity,  Eq.  (2.4)  can 
be  solved  by  a  perturbative  approach.  If  the  pulse  area  is 
increased  by  increasing  either  the  field  strength  or  the 
pulse's  duration,  the  perturbation  treatment  is  not  applic¬ 
able  and  there  is  no  analytical  solution  to  these  coupled 
differentio-integral  equations.  In  the  strong-field  limit 
(i.e  ,  y  »  T,y,  y  »ka,  and  \T  »  1 ,  where  T  is  a  pulse’s 
duration),  however,  it  may  be  possible  to  simplify  the 
QMTE  by  using  a  dressed-state  picture. 

B.  QMTE  in  dressed-state  picture 

Instead  of  using  the  bare-state  picture  above,  it  is  pos¬ 
sible  to  introduce  a  setniclassical  dressed-state  picture 
(DSP)  (Ref.  35)  in  which  the  dressed-state  basis 
represents  an  exact  solution  of  the  atom-field  interaction 
neglecting  relaxation.  The  relationship  between  dressed- 
state  density-matrix  elements,  denoted  by  pD  =  (pU0, 
Pbb’Pab'Pba^’  and  bare-state  density-matrix  elements, 
denoted  by  pg  =(pn,p22,pi2,p2l )  is  given  in  Appendix  A, 
along  with  the  corresponding  relaxation  parameters  in 
the  DSP.  The  time  development  in  the  DSP  is  given  by 

P,y{  v,f)  =  (L..  +  L  u  )pp{v,t ) 

h  f  w(  v'  —  v  )Kpn(v',t  i d  v'  ,  (2.6) 

where 


K  = 


cos  2(0‘  8)  sin:(0'—  6) 

sin;(0'  —  0)  cos2(0'  —  0) 

sin2  (0'  —  8)  —  sin2(0'  —  6) 


2  2 

sin2(fl'  —  0)  sin2(  0'  -  0) 
2  2 

1 


cos©=  — ^v7  1  +kv/2k  , 
V2 


sin@=  — 1  —kv/2X  , 
Vl 


sin2(  8'  —  0)  sin2  (O'  — 8) 

2  2 

sin2(t?'  —  0)  sin2(fl'  —  0) 
2  2 

cos2(0'  —  0)  stn2(0'—0) 
sin2(0'  —  0)  cos2(d'  —  0) 
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A  vo  1 '  ‘.  1  *  Lv  :< 
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and 


*=[r+(*v>  / 4] 


1/2 


The  matrix  LniJ  has  components 


1“*^.  33 

—  44 


where 

n/f=2A.  =  [4*2  +  (A:v)2]1/2 


(2.9) 


(2.10) 


(2.11) 


Matrix  elements  of  the  relaxation  matrix  Ly  are  given  in 
Appendix  A. 


III.  SOLUTION  OF  THE  QMTE  IN  DSP: 

STRONG-FIELD  REGIME 

The  dressed-state  bas:«  is  particularly  useful  '"her.  the 
generalized  Rabi  frequency  f lR  is  much  larger  than  the 
relevant  decay  parameters.  In  the  strong-field  limit,  an 
approximate  solution  to  Eq.  (2.6)  can  be  obtained  for  the 
difference  kernel  (2.5).  Assuming  that 

X»  all  Ps  ,  (3.1a) 

X»ka  (3.1b) 

(recall  that  a /V~l  is  the  rms  velocity  jump  per  collision), 
all  coupling  between  the  coherences  and  populations  in 
the  DSP  become  negligibly  small  (secular  approxima¬ 
tion36).  In  an  interaction  representation  defined  by 
Pab  ~Pabe  *  ■  the  dressed-state  density-matrix  elements 
_ I 


q(v,r)  = 


<7ao<v,r) 
<7w,<v,  t) 


1 

In 


fdv'I 


♦  -r.  fA(  v  •  v*)r  +  H  / 

are 

T. 


1  rip( 


( k  o  t  /  2  I* ) 


To  find  the  population  densities  paa  and  phb  from  q and 
qbh  using  Eq.  (3.5)  the  matrix  e  psSin  Eq.  (3.5)  needs  to 
be  calculated  analytically.  We  perform  this  calculation, 
assuming  that  the  system  is  closed  to  population  loss  via 
spontaneous  emission,  i.e.,  y,  =  0, y 2. i  ~Y2-  The  popula¬ 
tion  density  is  given  by 

ptla(v,n  =  Sllqga(v,n+Sl2qbh(v,n  ,  (3. Ha) 

pbblv,t  )=S  2iqaJv,t )+  Snqbh(v,t)  ,  (3.8b) 

where  the  are  given  by  Eq.  (B8)  in  Appendix  B. 

Dressed-state  coherence.  Before  solving  the  time- 
evolution  equation  quantitatively  for  the  coherence  densi¬ 
ty,  the  physical  properties  of  the  phase  shift  induced  by 
collisions  should  be  discussed.  The  signal  decay  resulting 
from  the  collision-induced  phase  shift  depends  critically 
on  the  phase  accumulated  in  the  time  interval  of  interest. 
This  phase  shift  is  measured  by  the  quantity 
i[n*<v')  — n„(v)]r|  appearing  in  Eq.  (3.2b).  If  the 
phase  shift  is  larger  than  unity,  the  contribution  to  the 


satisfy  the  uncoupled  time-evolution  equations 

Tj(v,t)  =  LpTj(v,/)+  Jdv'Tji  v',t )  ,  (3.2a) 

pabiv,t  )  =  l.yiipab(v,t)  +  f  dv'wiv'  —  v )pab ( v' , f ) 

X e'ln* * v  ’  1  v|r 


Pba  I  V,  f  )  ~  p*b  , 
where 


(3.2b) 

(3.2c) 


T)(  V.f  )  = 


Poa(v,t) 


and 


(3.3) 


(3.4) 


These  equations  may  now  be  solved  independently. 

Dressed-state  population.  First,  we  solve  Eq.  (3.2a)  by 
writing 

rf-e^'q  (3.5) 


and  substitute  it  into  (3.2a)  to  get 


q( v,r  )=  J d v'(e/p  '  ')  'eLp  q(  v',r )  . 


(3.6) 


Using  approximation  (3.1)  and  Fourier  transform  tech¬ 
niques,  one  can  obtain  the  integral  solution 


9aa* v  ’f0  ' 

Wv',6, ) 


(3.7) 


echo  s  .  •»  1  om  the  dipoles  is  reduced.  ThereTore,  it  is 
importa.i  r  nalyze  the  phase  shift  |[fl(v')  —  Il( v)]7!  in 
order  to  .  ihe  dominant  terms  and  to  neglect  terms  of 
order  of  o  /  \  or  a  /u.  We  expand 


(ll(v')-  ll(v)]r  = 


k  v  k  8vT 
[4v2  +  ( k  v)2]1 /2 
2\2(kbv)2T 
[4r+(*v)Ji,/J 


--  0(f 3  or  a /u  i 


(3  9) 


with  e  =  /cbv/y  and  bv  ~  |v~v''.  The  first  and  second 
terms  are  analyzed  for  three  velocity  ranges: 
k  v  «  \,k  v  =  \,k  v» In  the  analysis,  we  set  bv~-o. 

Consider  the  first  term  in  (3.9).  ( i )  If  kv<.<  y,  the  mag¬ 
nitude  of  the  first  term  is  approximately  k  vT  e.  By 
choosing  a  large  t  can  be  very  small  so  that  this  term 
vanishes,  (li)  If  kwy  or  kv»x .  the  magnitude  of  (iie 
first  term  is  of  order  kot.  As  the  pulse  duration  T  in¬ 
creases  such  that  koT  >  1,  no  matter  what  the  value  of 
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the  field  intensity  is,  the  phase  shift  is  significant  in  this 
time  interval.  Consequently,  collision-induced  decay  of 
atoms  having  kv>\  is  produced  for  time  intervals  T 
satisfying  kcrT>\,  even  with  large-field  intensities 
( T  <\<ku). 

The  magnitude  of  the  second  term  is  of  order  koT  e 
for  kv«\,  and  of  order  kbvT[  \(kcr  )2/(  ku  )3  j  for 
kv»x.  Hence,  for  sufficiently  small  e  and  y//ct/,  the 
second  term  in  (3.9)  is  smaller  than  the  first  one.  Howev¬ 
er,  this  does  not  necessarily  mean  that  the  contribution 
from  the  second  phase  change  term  can  be  neglected, 
since  the  absolute  value  of  the  second  term  may  become 
larger  than  unity  when  the  time  interval  is  extended.  For 
pulse  durations 

T  <  X-~  (3.10) 

( ka  T 

the  second  term  in  Eq.  (3.9)  can  be  neglected.  We  assume 
that  condition  (3.10)  is  satisfied  in  the  present  problem 
and  neglect  the  contribution  from  the  second  term  in  Eq. 

(3.9). 

Equation  (3.2b)  is  solved  with  the  phase  change 
[ftp(v')  — ft R(v)]r  replaced  by  the  first  term  of  the  ex¬ 
pansion  Eq.  (3.9).  Using  Fourier  transform  techniques 
and  the  strong-field  condition,  one  finds  formal  solutions 
for  the  coherence  densities 


IV.  CALCULATION  OF  EXTENDED-PULSE  PHOTON 
ECHOES:  STRONG-FIELD  REGIME 

The  dressed-state  solution  obtained  in  Sec.  Ill  is  used 
to  obtain  expressions  for  the  extended-pulse  echo  ampli¬ 
tude  for  the  following  three  cases:  (Sec.  IV  A)  Neglecting 
relaxation;  (Sec.  IV  Bi  including  relaxation  but  neglecting 
velocity-changing  collisions;  (Sec.  IV  C)  including  effects 
of  velocity-changing  collisions. 

A.  Neglecting  relaxation 

In  the  absence  of  relaxation,  the  strong-field  conditions 
(3.1)  are  always  satisfied.  In  this  limit,  the  photon-echo 
signal  can  be  calculated  by  standard  techniques.3'  The 
photon-echo  amplitude  depends  on  the  off-diagonal 
density-matrix  element  pi\(tc )  evaluated  at  the  echo  time 
te  +/,  —  /,  (see  Fig.  2).  For  an  initial  excitation  pulse 
of  short  duration  (kut]«\)  and  a  second  excitation 
pulse  of  arbitrary  duration  Tp2  =  ti  — 12  (see  Fig.  2),  one 
finds 

—  / 2  sin(  2 V i  T'i  )*2 

p2i(te) - -= - - 

V  TTU 

J  4^2  -F  ( Ac  v  )2 

(4.1) 


Pus(v,n~e 


f-vU  ‘  “lK  '[) 


X  f  dTdv'e 

J  X 

1  ~  ;!io  Ml  •  l„ 

Pha  (V.llsf  * 


ifclv-v’lW'lr.l-r  I  , 

Pablv’t 0>  - 
(3.11a) 


X 


f  d-d  v'e 


A  i  v  -  v  ir  4-  /  i  r.f  —  ) 

Pba^J  „)  , 


(3.11b) 


where 

/  (r,rl  =  T  f  dt  expj  ~-[(£t  ±r)ka/2]2\  (3.12) 

J  o 

and 


t--kv/ilR  . 

Equations  (3.8)  and  (3.1 1)  represent  solutions  to  QMTE 
in  the  DSP.  Several  approximations  were  made  in  arriv¬ 
ing  at  the  results;  in  particular,  it  was  assumed  that  (i) 
velocity-changing  collisions  act  in  a  state-independent 
manner;  (ii)  collisions  are  weak  and  characterized  by  a 
difference  kernel;  (iii)  spontaneous  emission  occurs  from 
level  2  to  level  I  only;  and  (iv)  the  Rabi  frequency  is 
greater  than  all  relaxation  parameters  and  is  smaller  than 
the  inhomogeneous  width  ( A: w )  of  the  atomic  velocity  dis¬ 
tribution. 

Solutions  (3.8)  and  (3.11)  can  be  transformed  back  to 
bare  state  to  get  corresponding  time-dependent  density- 
matrix  elements  that  are  needed  for  solving  transient 
problems  involving  excitation  pulses  of  long  duration. 


The  velocity  integration  is  computed  numerically.  The 
variation  of  the  echo  amplitude  |p2|U,  )|  with  duration  of 
the  second  pulse  is  plotted  in  Fig.  3  (the  corresponding 
result  for  the  weak-field  limit  is  shown  in  Fig.  10  of  Ref. 
31). 

In  contrast  to  the  weak-field  result  in  which  the  ampli¬ 
tude  increases  linearly  with  Tp2,  for  a  strong  excitation 
field,  the  echo  amplitude  oscillates  with  Tp2  as 

sin2|[^  +  (*v/2)J]l/Jrp2|  . 

In  the  strong-field  case,  the  velocity  range  of  atoms  con¬ 
tributing  to  the  echo  is  selected  by  the  Rabi  frequency  L2; 
that  is,  one  excites  atoms  having  kv<Xi ■  During  the 
second  pulse  interval,  the  phase  associated  with  atoms  in 
this  velocity  range  varies  approximately  as  X:Tp2.  Since 
*2^2  »1,  this  phase  is  significant  to  echo  formation.  It 
produces  an  oscillation  in  the  echo  amplitude  that  varies 
approximately  as  sin^T^.  The  signal  saturates  with  in¬ 
creasing  Tp2  as  the  population  difference  of  levels  1  and  2 
approaches  zero.  » 


0  t?  t3  t. 

FIG.  2.  Excitation  pulse  photon-echo  excitation  scheme. 
The  maximum  echo  signal  occurs  at  t,,  =  t ,  +  t:  - 1 , ,  denoted  by 
the  arrow. 
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In  the  weak-field  case,  since  Y2rp2«l,  the  velocity 
range  is  determined  by  the  condition  kv<(  \/Tpl)  rather 
than  kv<Xi-  In  a  first  approximation,  one  can  replace 


the  integrand  in  (4.1)  by  j(Tpl ),  while  restricting  the  ve¬ 
locity  range  to  ~(kTp2)~ As  a  result  the  echo  ampli¬ 
tude  increases  linearly  with  TP2- 


B.  Relaxation  of  echoes  neglecting  velocity-changing  collisions 


The  results  of  Sec.  IV  A  are  now  generalized  to  include  effects  of  spontaneous  decay  and  phase-interrupting  collisions. 
The  spontaneous  emission  rate  y2,  i  *s  taken  equal  to  y2  and  *s  taken  equal  to  y2/2  (“closed”  system  with  level  1  the 
ground  state).  From  time  t]  to  r2,  the  relevant  density-matrix  element  p,2  decays  exponentially  as 


p12(v,r2)  =  e 


=  „  IPph 


Pi2(v,t,  ) 


(4.2) 


In  the  second  field  region,  an  analytic  solution  is  no  longer  available  owing  to  the  decay  terms.  In  the  strong-field  limit, 
the  DSP  approach  discussed  in  Sec.  II  is  useful.  The  main  procedure  for  calculating  p2,(v,f3)  for  a  given  p12(v, t2)  in¬ 
volves  (a)  defining  initial  conditions  of  density-matrix  elements  in  the  DSP  [pD(v,t2 )]  in  terms  of  p,2(  v,r2 );  (b)  solving 
pD(v,t3 )  in  terms  of  pD(v,t2 );  and  (c)  transforming  pD(v,t} )  back  to  the  BSP  to  get  p2l(v,f3 )  using  Eq.  (Al)  of  Appen¬ 
dix  A.  When  this  procedure  is  carried  out,  one  may  obtain 


p21(v,/3)  = 


2^-'r<>1^ 
4y2  +  (kv)2 


exp 


trph-r2/2)Tp2 


(kv)2 


4y^  +  ikv)2 


-  cos)  [4y2  + (kv)2]  /2T  2 j  exp 


Ph  Yi/1)Tpi  2  2 


4y2  +  (k  vr 


p,2(v,t2)  , 


(4.3) 


where  p12(  v,f2 )  is  given  by  Eq.  (4.2).  In  the  interval  from 
t3  to  te  the  coherence  density  rephases  and  decays  at  the 
same  rate  as  that  in  the  second  interval.  Explicitly,  one 
finds 


By  combining  all  relevant  density-matrix  elements  to- 


FIG.  3.  Variation  of  the  photon-ccho  amplitude  with  the 
second  pulse's  duration  in  the  absence  of  relaxation  in  the 
strong-field  limit  In  this  as  in  Figs  4  and  5,  y/ku  =0.2  and 
A  «if.  —  1000. 


I 

gether  and  integrating  over  velocity  numerically  at 
t  =/„  =z3  +  t2  —  1 1 ,  one  is  able  to  depict  the  variation  of 
echo  amplitude  with  the  second  pulse’s  duration  in  the 
presence  of  relaxation  (see  Fig.  4).  It  is  seen  that  the  echo 
signal  oscillates  and  increases  with  V  In  this  case, 
there  are  two  factors  responsible  for  the  increase  of  echo 
amplitude.  The  first  factor  has  been  already  discussed 
above  for  the  case  without  relaxation;  namely,  the  in- 


FIG.  4.  Variation  of  the  photon-echo  amplitude  with  the 
second  pulse's  duration  in  the  absence  of  velocity  changing  col¬ 
lisions  in  the  strong-field  limit  (rph/y  =  0. 1  and  y;/*  =  0.01 ). 
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crease  of  the  pulse  area  excites  more  dipoles  to  emit  be¬ 
fore  saturation  is  reached.  The  second  factor  is  related  to 
the  fact  that,  in  each  time  interval,  the  density-matrix  ele¬ 
ment  contributing  to  echo  formation  decays  at  a  rate  that 
depends  on  whether  the  density-matrix  element  is  associ¬ 
ated  with  level  population  or  coherence.  In  the  two 
field-free  regions,  the  coherence  p12  or  p2,  decays  at  rate 
rph  +  yl2.  In  the  second  field  region,  the  second  pulse  of 
strong  intensity  continually  switches  atoms  between  pop¬ 
ulation  and  coherence,  and  the  decay  rate  in  this  region 
can  be  associated  with  coherence  (at  rate  rph  +  y12)  or 
population  (at  rate  y2).  As  the  pulse  duration  is  in¬ 
creased,  the  population  decay  is  increased  relative  to 
coherence  decay.  Consequently,  if  rph  +  y  ,2  >  y2,  the  in¬ 
creased  population  time  implies  a  decrease  in  the  echo 
decay  rate  with  increasing  duration  of  the  second  pulse 
(recall  that  a  decrease  of  decay  rate  implies  an  increase  in 
echo  amplitude)  when  increasing  the  second  pulse’s  dura¬ 
tion.  If  rph  +  y|2<y2  (only  possible  at  low  perturber 
pressures),  the  situation  is  reversed.  This  feature  das 
been  also  noted  in  the  weak-field  result. 

It  is  interesting  to  compare  schematically  the  echo  de¬ 
cay  rate  (to  be  defined  more  precisely  below)  predicted  by 
Eqs.  (4.3)  and  (4.4)  in  Fig.  5  with  the  experimental  one  of 
Yodh  et  al."  in  Fig.  6,  even  though  the  result  does  not 
yet  include  effects  of  velocity-changing  collisions.  In 
comparing  Figs.  5  and  6,  one  sees  that  the  phenomenon 
common  to  both  is  that  the  decay  rate  decreases  with  in¬ 
creasing  the  second  pulse  interval.  The  results  differ  in 
that,  in  Fig.  5,  the  decay  rate  decreases  at  the  same  rate 


FIG  5.  Considering  spontaneous  emission  and  phase¬ 
interrupting  collisions  only,  one  finds  that  the  echo  decay  rate 
envelope  decreases  at  a  constant  rate  for  al!  Tpl,  which  differs 
from  the  behavior  of  the  echo  decay  rate  in  the  presence  of 
velocity-changing  collisions.  (An  oscillation  of  f}  with  Tpi  is 
seen  in  this  result,  which  has  not  been  averaged  over  a  distribu¬ 
tion  of  field  strengths. )  In  this  figure,  rph/*=0.01  and 
Yi/X^O.Ol. 


for  all  values  of  Tpl,  whereas,  in  the  experiment,  the  de¬ 
cay  rate  decreases  very  slowly  for  small  T,2  but  more  rap¬ 
idly  for  Tp,  comparable  to  the  observation  time  of  the 
echo.  This  indicates  that  spontaneous  emission  and 
phase-interrupting  collisions  alone  do  not  provide  an  ap¬ 
propriate  physical  picture  for  interpreting  the  experimen¬ 
tal  result. 

C.  Effects  of  velocity-changing  collisions  on  photon  echoes 

In  order  to  present  a  theory  of  extended-pulse  photon 
echoes  including  effects  of  velocity-changing  collisions 
which  is  relevant  to  the  experiment  of  Yodh  et  al.,  we 
first  review  some  of  the  parameters  of  that  experiment. 
The  active  atom  is  l74Yb,  the  transition  is 
[(6s2)'S0  —  (6.s6p  i’F,  ]  having  wavelength  X  =  555.6  nm, 
and  the  perturber  atoms  are  argon.  The  echo  time 
/,  =  1200  nsec  is  fixed  throughout  the  experiment.  The 
second  pulse’s  duration  Tp2  varies  from  40  —  960  nsec. 
For  each  value  of  TP»  the  echo  intensity  is  measured  as  a 
function  of  argon  pressure,  varying  from  10-40  mtorr. 
The  decay  rate  is  calculated  from  each  measurement  and 
is  plotted  as  a  function  of  Tp2,  with  the  expeiimental  re¬ 
sults  indicated  by  crosses  in  Fig.  6.  When  Tp2  is  less  than 
480  nsec,  the  curve  decreases  slowly.  When  Tp2  is  close 
to  te ,  the  curve  decreases  more  rapidly. 

The  temperature  of  Yb  vapor  is  constant  at  500  °C. 
The  first  pulse  area  (as  defined  by  2y,  T  , )  is  of  order  uni¬ 
ty.  The  Rabi  frequency  of  the  second  excitation  field  is 


FIG.  6.  Variation  of  the  decay  rate  (/?)  of  the  photon  echo 
with  the  second  pulse's  duration  Tp2.  The  curve  represents  the 
theoretical  result  for  the  strong-field  limit  and  the  crosses 
represent  the  experimental  data.  They  are  normalized  at 
Tpi/t, -0.23.  A  Gaussian  distribution  for  y  with 
Xm,,/ku  =0.28  has  been  incorporated  into  the  theory.  The  oth¬ 
er  parameters  are  chosen  to  correspond  to  the  experimental 
ones  given  this  value  of  yma,,  namely,  r /v-«,'5x  io  ’. 
y/X -  !  ■  1  x  1 0  \  and  k  <  o  /ym„  1  =  64X10  2 
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y~10Vs  (in  the  following  discussion,  y  indicates  the 
Rabi  frequency  of  the  second  field).  The  natural  lifetime 
of  the  upper  level  is  875  nsec  and  the  corresponding 
spontaneous  emission  rate  is  y2  =  1. 14  X  lOVs.  The  total 
cross  sections  for  phase-interrupting  collisions  and 
velocity-changing  collisions  avcc  are  306  and  396  A  \  re¬ 
spectively.  The  velocity-changing  collision  kernel  width 
is  6v~0. 57  m/s  [(kdv=(2n/k)u -6.4X  lOVs],  The 
relative  average  speed  is  v,  ~  1000  m/s.  The  average 
speed  of  Yb  is  v  —  306  m/s.  The  most  probable  speed  of 
Yb  is  u  —  270  m/s  f  ku  =  ( 2n/k)-~  3  X  10'‘/s].  The  re¬ 
duced  mass  is  p  ~  27.  24  X  10  27  kg;  the  velocity-changing 
collision  and  phase-interrupting  collision  rates  are 
r  =  ,Vv,CTM.(  -0.5X  106/s  (for  10  mtorr)  and  Tph 
-  jVvroph  —  0. 39  X  lOVs  (for  10  mtorn,  respectively.  It  is 
necessary  to  examine  some  ratios  of  these  parameters, 
since  they  are  relevant  to  the  approximations  in  the  cal¬ 
culation.  The  strong-field  limit  requires  that  the  Rabi 
frequency  be  greater  than  all  relaxation  rates,  which  is 
satisfied  here,  i.e.,  T/y~5X10  ’,  yph/y -3.7X10 
y/y~  1.1X10  ',  ko/y~-  6.4X  10  2.  We  keep  terms  to 
first  order  in  these  ratios  in  the  following  calculations. 
Collisional  relaxation  rates  T  and  Fph  and  spontaneous 
emission  rate  y2  are  all  of  the  same  order  of  magnitude; 
hence  all  should  be  considered  in  the  calculation.  Quan- 
_ I 


tities  T  Tpl  and  koTpl  are  of  the  order  of  unity  for  large 
Tpl  so  that  effects  of  velocity-changing  collisions  cannot 
be  neglected  in  the  second  pulse  interval,  since  7.,  varies 
over  a  large  range.  The  pulse  area  \Tpl  is  larger  than 
unity  for  all  values  of  Tp!  in  the  experiment  so  that  a  per¬ 
turbation  approach  fails.  The  ratio  of  Rabi  frequency  to 
inhomogeneous  width  ku  is  on  the  order  of  0. 1 .  Terms  of 
second  or  higher  order  in  y /ku  will  be  neglected.  The  in¬ 
equality  kuTpl'>>\  holds  for  all  values  of  Tp,  in  the  ex¬ 
periment,  a  condition  that  simplifies  the  integration  in¬ 
volved  in  the  Fourier  transformation  in  the  second  pulse 
interval. 

For  the  above  experimental  parameters,  we  evaluate 
the  extended-pulse  photon-echo  amplitude  in  the  pres 
ence  of  spontaneous  emission,  phase-interrupting  col¬ 
lisions  and  velocity-changing  collisions.  The  density- 
matrix  elements  contributing  to  echo  formation  are  cal¬ 
culated  in  each  time  interval. 

In  the  first  time  interval  (0  » / ,  ),  the  pulse  area  is 
larger  than  unity  and  a  perturbation  solution  is  not  valid. 
Since  T Tp ,  (or  yTpli  «  1 ,  and  kaTP\  «  I,  the  relaxation 
in  this  region  can  still  be  neglected.  In  the  second  inter¬ 
val  (/,-/,),  there  is  no  excitation  field  and  the  coherence 
density  is  the  same  as  that  derived  in  the  weak-field  limit. 


Pl2(v2-f2  >=  ~ 


.  _  _  if  ♦  f'  K  /Av,  »i/, 

i  sin(2ylT„|  )e  r  •  ‘ 


4ttV  7 TU 


f  /  ‘  ' dr dv e‘k  ':  "r 
I  r'-’  b 

X  exp  T  J  dl 'expj  -  [kai  t '  +  r)/2]‘  | 


(4.5) 


The  present  task  is  to  solve  for  p2l(\’,l,)  in  terms  of  p } ,( v,  f  ? )  in  the  third  time  interval  ((,  /,),  since  it  is  this  term 

which  gives  rise  to  the  echo  signal  having  a  maximum  amplitude  at  i  I,.. 

One  can  use  the  strong-field  solution  of  the  dressed  state  OM  I  F  obtained  in  Sec.  Ill,  and  transform  it  back  to  the 
bare-state  representation  to  obtain  p2)(v,r,).  The  process  for  doing  this  has  already  been  carried  out  in  part  B  of  this 
section.  By  modifying  the  calculation  to  include  the  effects  of  velocity-changing  collisions,  one  arrives  at 


P 2l<vl.'l  >  = 


sin20, 


SM  Si,  5p  IS,,)  f  </v,£/rsin2W,e 
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A  <  V  ,  %  ,  1  r  *  /  '  t  r,  / 


(  cos' 


1 

2  rr 


J  dv2d  r  sin' (be 


<r  •  i  '  -.  I 


(4.6) 


where  the  Vs  are  given  in  Appendix  B  and  /  is  defined  in  Fq  (3.12).  Note  that  f>,  is  a  function  of  \\ . 
In  the  last  interval  ((,  f, ),  the  calculation  is  similar  to  that  for  p  |2(v,(, ),  one  finds 


Pi  f  ^  ^  ^ 


'f  *  l  ph  *  (AV4.U  t  y 


2tt 


//"  d  t  d  v,e  '  *  exp  I  V  J"  dl '  exp  j  [  k o(  f '  r  1/2  ]‘  |  L»,,(  v,,  t ,  1 


(4.7) 


By  carrying  out  the  velocity  integrations  at  t  +r2  —  r , ,  we  find 
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^  1  i  r  *-  r  h  f  y,,H +  <2TV  rr/ko  )crfl  k<7t,  /2>  r  +  ■* 

'’""•’“'THT 


dv 


f  — i  ~  — , !  CDP  +  CDC  )  , 

J  ■«  4*2  +  (*v)2  UP  Dt 


CDP=  exp 


rr,2exp[-(*af2/2)2]+(rph-y2/2>r,2 


CDC=  -  cosj(4Af  +(kvr]  /-7p2|exp 


K2//2)7'p2  , 


_ I  fc  vF 

4v2  +  <^v)2 

2l!__ 


4.8) 


4;y~  +  i  k  v  i 


+  exp(  W  -rr/tjka  )  |  e  r  f[  <  £:  7^  2  + r  2  )&  tr/2]  —  erf(  /c  err  2  /2 )  | 
+  exp(  TV v/cka  ) j er{[(^Tp2~ t2  )kcr  /2]~  erf(kot2/2)\  , 


where  5  is  defined  in  Eq.  (3.12).  The  symbols  CDP  and 
CDC  denote  contributions  from  dressed-state  population 
and  coherence,  respectively. 

This  is  an  approximate  analytic  solution  for  an 
extended-pulse  photon  echo  assuming  strong  excitation 
fields.  In  order  to  compare  this  result  with  the  experi¬ 
ment,  we  define  'he  decay  rate  of  the  echo.  The  echo  in¬ 
tensity  is 


,  ,  |2  -»Vi> 

h-h)Pi\\  -e 


(4.9) 


where  P  is  the  perturber  pressure  and  /?  is  referred  to  as  a 
decay  rate  which  is  a  function  of  Tp1’  defined  by 


/?( Tp2 ) 


log/,//,, 

P 


(even  though  the  ft  does  not  have  dimensions  of  a  rate,  it 
is  defined  as  a  decay  rate  to  be  consistent  with  the  nomen¬ 
clature  of  Ref.  1 1).  The  decay  rate  j3  versus  T,i  for  the 
strong-field  limit  is  plotted  in  Fig.  6.  A  Gaussian  distri¬ 
bution  of  the  Rabi  frequencies  (to  reflect  the  spatial  prop¬ 
el  ,iu  uf  the  light  beam)  his  L.::r.  included  in  '.he  numeri¬ 
cal  result,  with  a  maximum  value  of  taken  equal 

to  0.28.  Our  theoretical  result  (represented  by  the  curve) 
fits  the  experimental  data  (represented  by  the  crosses) 
very  well. 

It  should  be  noted  that,  although  Eqs.  (4.8)  and  (4.9) 
have  been  compared  with  experiment,  they  do  not  corre¬ 
spond  exactly  to  the  quantity  measured  in  the  experi¬ 
ment.1  The  signal  at  t—t,  actually  consists  of  three 
parts.'11  First,  there  is  the  “true”  echo  contribution 
which  peaks  at  t  =  te  and  arises  from  the  dependence  of 
I  F2|(G(  on  Pt2*G).  This  is  the  term  calculated  in  this  pa- 

I  per.  Second,  there  is  a  free-induetton  decay  (FID)  signal 

emitted  after  the  second  pulse  which  arises  from  the 
!  dependence  of  p2,(  r, )  on  p,,(  t2 1  and  p22l  >2  )•  Third,  there 

I  is  an  FID-like  signal  emitted  after  the  second  pulse  which 

k  arises  from  the  dependence  of  p2 1^3)  on 

experiment,  a  random  phase  was  introduced  between  the 
two-excitation  fields  and  the  signal  was  measured  as  the 
>  difference  in  signals  produced  with  and  without  the  first 

field  present.  This  procedure  was  intended  to  eliminate 
the  FID  background.  It  can  be  shown  that  such  a  pro¬ 
cedure  does  not  completely  eliminate  the  contributions  to 
the  signal  from  the  second  and  third  terms  discussed 
above.  Under  those  experimental  conditions,  however,  it 


can  be  estimated  the  reported  signal  is  approximately 
represented  by  the  echo  contribution  [4.8)  and  (4.9)]. 

To  further  pursue  some  of  the  implications  of  the  re¬ 
sults,  it  is  useful  to  compare  the  strong-  and  weak-field 
solutions,  and  to  determine  which  components  of  the 
solution  contribute  to  the  dependence  of  the  echo  decay 
rate  on  Tp2. 


V.  ANALYSIS  OF  THE  RESULTS 

From  the  strong-field  result,  we  find  that  the  decay  rate 
of  the  echo  amplitude  decreases  with  increasing  the 
second  pulse’s  duration,  which  is  in  good  agreement  with 
the  experimental  results  of  Yodh  el  al."  Although  we 
have  considered  a  strong-field  limit,  the  qualitative  ex¬ 
planation  of  the  dependence  of  P  on  Tp2  remains  the  same 
as  in  the  weak-field,  perturbative  domain.  In  discussing 
the  result,  it  is  always  assumed  that  kate>\  and 
kuTpl  >  1. 

We  first  consider  the  limit  that  ka  Tp2<\,  7  Tp,  >  1 . 
For  koTpl<  1,  the  time  interval  between  pulses  1  and  2 
and  between  pulse  2  and  the  echo  is  such  that 
kaitt,  —  Tpl)  ">  1.  This  collision-induced  phase  shift  is 
sufficiently  large  to  destroy  the  contribution  to  the  echo 
signal  from  any  atom  that  has  undergone  a  collision. 
Consequently  the  echo  intensity  decays  homogeneously  as 

exp[  —  (y ,,  +  rph  +  T  )te  +  2W  Tt/ka  eT([kate  /4 )  ]  , 

independent  of  Tp2 ■  There  is,  however,  a  weak  depen¬ 
dence  of  the  log  of  the  echo  intensity  on  V  As  was 
shown  in  the  weak-field  limit,51  an  increase  in  TP2  results 
in  an  increase  in  the  number  of  atomic  dipoles  contribut¬ 
ing  to  the  echo  signal,  provided  that  kaTp2<\.  If 
kuTp2>  1,  the  echo  intensity  varies  as  (  X2Tp2r,  implying 
that  log(///0)  varies  as  2  log(  7~p2 ).  Thus,  /3(  Tpl )  as 
defined  in  Eq.  (4.9)  decreases  slowly  with  increasing  Tpl 
when  ka  Tp2<  1.  For  intense  fields  such  that  X:Tp 2>  1, 
this  dependence  is  modified  somewhat.  On  averaging 
over  a  distribution  of  field  strength,  one  still  finds  a  slow 
decrease  in  p(Tpl)  with  increasing  Tp2  when  koTp ;<1 
and  yTpl<  1. 

As  Tp2  increases,  the  time  interval  ( —  Tp2)  decreases. 
For  sufficiently  large  Tp2,  one  arrives  at  the  limit  that 
ka(tc  T  2)«\,  such  that  velocity-changing  collisions 
result  in  negligible  phase  shifts  between  the  pulses.  Th^ 
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role  of  collisions  during  the  second  pulse  now  becomes 
important  since  ka  T,2>  1.  As  has  been  discussed  in  the 
weak-field  case,21  collisions  tend  to  limit  the  contribution 
irom  the  second  pulse;  only  the  beginning  and  end  of  the 
pulse  interval  contribute  significantly  to  echo  formation. 
In  effect,  the  EPPE  becomes  equivalent  to  a  series  of 
stimulated  photon  echoes  which  are  selected  in  a  manner 
that  renders  the  role  of  velocity-changing  collisions  negli¬ 
gible.  Thus,  during  part  of  the  second  pulse  (of  approxi¬ 
mate  duration  V-  there  is  no  collisional  decay  of  the 


signal.  Consequently^  the  signal  decays  as 
exp[  —  T(f,  —  Tpl)]  =  e  'e  p 2  and  log(///0)  varies  as 
T Tpl.  Thus,  0  varies  as  —  YTpl/P  for  large  Tp2.  [There 
is  an  additional  amplitude  factor  (T7'p,)  in  the  intensity 
which  slightly  modifies  this  result.]  In  going  from  weak 
to  strong  fields,  one  finds  that  the  amplitude  of  the  echo 
intensity  is  modified  somewhat,  but  that  the  overall  quali¬ 
tative  dependence  of  log///0  on  Tp2  is, fairly  insensitive 
to  the  field  strength.  Thus,  as  Tp2  is-increased  0(Yp2) 
goes  from  a  —  loglTT^)  to  a  —TTpl  dependence.  (An 
additional  dependence  of  0  on  Tp2  arises  from  the  fact 
that  population  and  coherence  decay  at  different  rates. 
This  results  in  a  factor  exp(2(  Tph  —  y2)Tp2]  in  the  echo 
intensity  which  also  contributes  to  0. ) 

Our  explanation  of  the  dependence  of  0  on  Tp2  differs 
from  that  of  Yodh  et  a/.,11  who  based  their  explanation 
on  a  strong-field  quenching  of  collisional  effects.  The  ar¬ 
guments  given  by  Yodh  et  al.  imply  that  strong-field 
quenching  becomes  important  when 


Tp2<  1 


a  condition  reached  in  their  experiment.  In  examining 
expansion  (3.9),  however,  we  find  that  the  appropriate 
condition  for  strong-field  quenching  is 


of  photon-echo  formation  and  relaxation  effects  arising 
from  spontaneous  emission,  phase-interrupting  collisions, 
and  velocity-changing  collisions.  The  dressed-state- 
picture  approach  has  been  used  for  obtaining  an  analyti¬ 
cal  expression  for  the  echo  amplitude  as  a  function  of  the 
second  pulse’s  duration  in  the  strong-field  limit.  The 
techniques  employed  in  this  paper  may  be  used  to  study 
the  collisional  modification  of  intense  radiation-matter  in¬ 
teractions  in  coherent  transient  spectroscopy.  We  have 
also  developed  numerical  methods  for  solving  the  QMTE 
for  arbitrary  field  intensity.40  The  numerical  results  are 
in  good  agreement  with  the  analytical  ones  in  both  the 
weak-  and  strong-field  limits.  Moreover,  such  methods 
can  be  further  applied  to  problems  involving  more  realis¬ 
tic  collision  kernels. 
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APPENDIX  A 

The  transformation  of  density-matrix  elements  ex¬ 
pressed  in  bare-state  and  dressed-state  pictures  is  given 
below, 

Pu<,  =  T[pn+p2’  +  (Pii-p>2>cos20+(Pi2+P2i>sin20]  , 


Pbb  =  y[pn  +  P22 "<Pi i  Pit >  cos2P  —  (p,2+Pji  >  sin20]  , 


PU*  =  T[(P22-Pn>s'n20+P!2”P2l+<Pl2  +  P21>^s20]  - 


[4^2  +  (kv)2]l/2 

For  ka  TP>>  1  and  x  «ku,  inequality  (4.11)  is  never 
achieved  for  atoms  having  velocities  kv>\  for  which 
collision-induced  decay  cannot  be  quenched.  In  effect, 
Yodh  et  al.  implicitly  used  a  simplified  theoretical  model 
in  which  the  first  excitation  pulse  was  sufficiently  weak 
and  long  to  excite  a  narrow  range  of  atomic  velocities 
kv~{Tn)  «X2-  For  such  an  excitation  scheme,  Eq. 
(4.10)  would  be  the  correct  condition  for  strong-field 
quenching;  however,  in  the  actual  experiment,  the  first 
excitation  pulse  was  of  sufficient  intensity  (X\~ X2/l)  to 
require  the  more  detailed  analysis  given  in  this  work. 

VI.  CONCLUSION 

This  paper  aims  to  present  a  systematic  theoretical 
analysis  of  the  phenomenon  of  photon  echoes.  The 
theoretical  results  are  shown  to  be  in  good  agreement 
with  the  photon-echo  experiment  carried  out  by  Yodh 
et  al . "  A  physical  interpretation  of  the  experimental  re¬ 
sult,  that  the  photon-echo  decay  rate  decreases 
significantly  with  the  second  pulse’s  duration,  was 
developed,  based  on  the  study  of  the  fundamental  theory 


Pba  =  y[<P22~Pu'sin20-p,2-p,,  +  <Pi2+P2i  )cos2G]  , 


Pi  i  -  HP*  +P/>fr  +(P™  Pbb  >  cos2 0~(pul,  +pha  )  sin20]  , 


Pn-\[pua+Phh-(Puo  ~Pbh)cos20  +  {patl+PfH1  )sin20]  , 


Pi2=  \\lPuu  ~Pbb  >  sin2 0+puh -Pha  +(pah  +  Phu  )  cos20]  , 


P21  =  y [<P00  Pbb )  sin20+p6u  -poh  +  ( pab  +Pha )  cos20]  , 


where  sinP  and  cos 6  are  defined  in  Eqs.  (2.8a)  and  (2.8b). 

The  matrix  elements  of  the  relaxation  term  appearing 
in  Eq.  (2.6)  of  the  paper  are  given  below, 

r  _  sin220  ^  4n  „  .  4„ 

Ey| | -y2J — - - Tn  cos  0—  r,,stn  0 


■(r.j  +  r. 


sin220 


6308 


RU-WANG  SUNG  AND  PAUL  R.  BERMAN 


39 


■^*y!2  ^2,1  COS40-(r  1 1  +  f22  ) 


sin220 


+<r12+r21) 


sin22g 
4  ’ 


(A2b) 

iri3  =  t(>,2.i  +  rti_ri2)cos2<9  +  (r21-r22)sin253^y^  , 

(A2c) 

^yu=[(y2.1+r„-r21)cos20+(r12-r22)sin20]i^> 

<A2d) 

^y2.=y2.isin46»-(r11  +  r22)^^+(rI2+r21)sin22(9 , 


L  y22  =  y  2, ,  -  r , ,  sin40  -  r  22  cos *e 

sin220 


-<rl2+r2I) 


4 


L  2  24  ”  [  ( y  2>  1 4* r ||  T i2) sin20 4- (  T2,  —  T22 ) cos20] 
i  =  [( r  n  —  r  12)  cos20  — (y2I  +  T22  —  T21 )  sin20] 


^v34=(r12+r21-rll-r22-72.1)^^  , 

^y4i  =  [(r,1-r2l)cos20-(y2J  +  r22-r2,)sin20) 


(A2k) 


(A21) 


2  ,  sin20 


2 

(A2m) 


Ty42— [^11  r,2)sin20  ( y  2,  l  T22  r21)COS0]  ^  > 

(A2n) 

sin220 


Ly43 — ( r12+ r2i  I'd  r22  y2i ])■ 


4  ’ 


(A2o) 


lyAA=  -(rn+r22+y2il)^~-r2lcos*e-rllSm*e , 


rn  =  r  +  yi . 
r22=r +y2) 
r,2=r2,  =  r+rph+y12 . 


(A3a) 

(A3b) 

iA3c) 


APPENDIX  B 


The  evaluation  of  the  exponential  matrix  ea'=S  is 
given  in  this  appendix.  The  matrix  G 


G  = 


*11  ~8  22 
—  *11  8  22 

has  elements 

sin220 


(A2e) 

(A2f) 

^r23=[(y2.i  +  r'll-r21)sin20+(r12-r22)dos2i9]^— • , 

(A2g) 

,  sin20 
2  ’ 

(A2h) 

sin20 
2  ’ 

(A2i) 

^v32=[(rn-r2I)sin20-(y2  ,  +  r22-rl2)cos20]-51—  , 

(A2j) 

Lr»  =  -< rn  +  r22+y,,)^ - r12 cos4e- r21  sin4@  , 


*11“ - 
822 =  rp 


p  2 

sin22f? 


+  y2  sin40 


4-y2cos40  , 


(Bl) 

(B2a) 

(B2b) 


for  y,=0  and  y2  ,  =  y2.  To  diagonalize  the  matrix  G,  we 
find  the  eigenvalues  of  G  as  2.,=0  and  A2=gu  4-g22,  and 
define  a  matrix 


U- 


“11  “21 


[“  12  “22 
with  elements 
u 


8  22 


"  l82n+8222)'/2  ' 

*n 

<g?.+Si2>1/J  ’ 

Uli=~^2  =_“22  ‘ 
The  inverse  matrix  t/_l  is 

f  “  21  “21 


U~  —D 
where 


“12  _“l! 


v/2(g2,+g22)l/2 

*11  +*22 
Using  the  transformation 


D  = 


0  /2 


U 


- 1 


the  four  elements  of  S  can  be  found  as 


Sn  e 


_ -n 


*22 


*1!  A,/ 


_ -n 


S|2  e 


S2,  —e 


*II~*'*22  *ii"T*22 


g,!  (i -«*■■) 


-n 


*11  '*■*22 
*11 


*11  +*22 


(A2p) 

S22=e~r' 

*11 

1  822 

where 

*1!  +*22 

*11  +*22 

(B3) 

(B4a) 

(B4b) 

(B4c) 

(B5) 

(B6) 

(B7) 

(B8a) 

(B8b) 

(B8c) 

(B8d) 


39 


THEORY  OF  OPTICAL  COHERENT  TRANSIENTS  INCLUDING  .  .  . 


6309 


’Present  address:  Department  of  Physics,  Beloit  College,  Beloit, 
WI  5351  1. 

’W.  Demtroder,  Laser  Spectroscopy  (Springer,  Berlin,  1981). 

2M.  D.  Levenson,  Introduction  to  Nonlinear  Laser  Spectroscopy 
(Academic.  New  York,  1982). 

7R.  Shoemaker,  in  Laser  and  Coherence  Spectroscopy,  edited  by 
J.  T.  Steinfeld  (Plenum,  New  York,  1978),  p.  197. 

4I.  D.  Abelta,  N.  A.  Kurnit,  and  S.  R.  Hartmann,  Phys,  Rev. 
141,391  (1966). 

5M.  Scully,  M.  J.  Stephen,  and  D.  C.  Burnham,  Phys.  Rev.  171, 
213  (1968). 

6R.  E.  Slusher  and  W.  J.  Tomlinson,  Phys.  Rev.  179,  294  (1969). 
7J.  Schmidt,  P.  R.  Berman,  and  R.  G.  Brewer,  Phys.  Rev.  Lett. 
31,1103(1973). 

8P.  R.  Berman,  J.  M.  Levy,  and  R.  G.  Brewer,  Phys.  Rev.  A  11, 
1668(1975). 

9P.  R.  Berman,  New  Trends  in  Atomic  Physics  (North-Holland, 
Amsterdam,  1984),  p.  453. 

I0P.  R.  Berman,  Advances  in  Atomic  and  Molecular  Physics 
(Academic,  New  York,  1977),  Vol.  13,  p.  57. 
llA.  G.  Yodh,  J.  Golub,  and  T.  W.  Mossberg,  Phys.  Rev.  Lett. 
53,659(1984). 

12M.  Yamanoi  and  J.  H.  F.berly,  J.  Phys.  Soc.  Am.  B  1,  751 
(1984). 

13E.  Hanamura,  5.  Phys.  Soc.  Jpn.  52,  2258  (1983);  52,  3678 
(1983). 

,4N.  G.  Van  Kampen,  Stochastic  Processes  in  Physics  and  Chem¬ 
istry  (North-Holland,  Amsterdam,  1981). 

15P.  R.  Berman  and  W.  E  Lamb,  Phys.  Rev.  A  2,  2435  (1970). 
I6P.  R.  Berman,  and  W.  E.  Lamb,  Phys.  Rev.  A  4,  319  (1971). 
17P.  R.  Berman,  Phys.  Rev.  A  5,  927  (1972).  - 
18P.  R.  Berman,  Phys.  Rep.  43,  101  (1978). 

19P.  R.  Berman,  J.  Opt.  Soc.  Am.  B  3,  572  (1986). 

20P.  R.  Berman,  J.  Opt.  Soc.  Am.  B  3,  564  (1986). 

21M.  S.  Feld  and  A.  Javan,  Phys.  Rev.  177,  540  1 1969). 

22P.  R.  Berman,  P.  F.  Liao,  and  J.  E.  Bjorkholm,  Phys.  Rev.  A 


20,  2389(1979). 

21R.  Salomaa  and  S.  Stenholm,  J.  Phys.  B  8,  1795  (1975);  9,  1221 
(1976). 

24R.  M.  Whitley  and  C.  R.  Stroud,  Phys.  Rev.  A  14,  1498 
(1976). 

2SC.  Cohen-Tannoudji  and  S.  Reynaud,  J.  Phys.  B  10,  2311 
(1977). 

2fiP.  R.  Berman  and  R.  Salomaa,  Phys.  Rev.  A  25. 2667  (1982). 

27J.  E.  Thomas,  R.  A.  Forber,  L.  A.  Spinelli,  and  M.  S.  Feld, 
Fhys.  Rev.  Lett.  51,  2194  (1983). 

28T.  W.  Mossberg,  R.  Kachru,  and  S.  R.  Hartmann,  Phys.  Rev. 
Lett.  44,  73(1980). 

29 A.  Flusberg,  Opt.  Commun.  29,  123  (1979). 

20R.  A.  Forber,  L.  Spinelli,  J.  E.  Thomas,  and  M.  S.  Feld,  Phys. 
Rev.  Lett.  50,  331  (1983). 

31 R.  Sung  and  P.  R.  Berman,  preceding  paper,  Phys.  Rev.  A  39, 
6284(1989). 

32S.  Kryszewski  and  G.  Nienhuis,  J.  Phys.  B  20,  3027  (1987). 

3JJ.  Keilson  and  K.  E.  Storer,  Q.  Appl.  Math.  10,  243  (1952). 

34R.  F.  Snider,  Phys.  Rev.  A  33,  178  (1986). 

35E.  Courtens  and  A.  Szoke,  Phys.  Rev.  A  15,  1588  (1977). 

36C.  Cohen-Tannoudji,  in  Frontiers  in  Laser  Spectroscopy,  edited 
by  R.  Balian,  S.  Haroche,  and  S.  Liberman  (North-Holland, 
Amsterdam,  1977),  Vol.  1,  p.  3. 

37See,  for  example,  T.  W.  Mossberg,  R.  Kachra,  and  S.  R. 
Hartmann,  Phys.  Rev.  A  20,  1976  (1979). 

3xIn  the  experiment  of  Yodh  et  al.,  the  inequality  kuTpl  «  1  is 
not  satisfied.  The  dependence  of  echo  amplitude  on  the  first 
pulse  area  is  neglected  in  this  paper.  The  fact  that  kuTpX  >  1 
implies  that  the  first  pulse  excites  only  those  atoms  with  ve¬ 
locities  satisfying  kv<x  (assuming  X?pX  ~  1 ).  This  would  not 
significantly  modify  the  results  since  the  second  pulse  is  al¬ 
ready  providing  a  similar  velocity  selectivity. 

•WN.  Lu  (unpublished). 

■^R.  Sung,  Ph.D.  thesis  (New  York  University,  University 
Microfilms,  1988)  (unpublished). 


